Non-contact gears: II. Casimir torque between concentric corrugated cylinders for the 

scalar case 
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The Casimir interaction between two concentric corrugated cylinders provides the mechanism 
for non-contact gears. To this end, we calculate the Casimir torque between two such cylinders, 
described by (5-potentials, which interact through a scalar field. We derive analytic expressions for 
the Casimir torque for the case when the corrugation amplitudes are small in comparison to the 
corrugation wavelengths. We derive explicit results for the Dirichlet case, and exact results for the 
weak coupling limit, in the leading order. The results for the corrugated cylinders approach the 
corresponding expressions for the case of corrugated parallel plates in the limit of large radii of 
cylinders (relative to the difference in their radii) while keeping the corrugation wavelength fixed. 



I. INTRODUCTION 



C^' The Casimir torque between two material bodies, which is the rotational analog of the Casimir force was studied 
^ for the first time in 1973 0. The Casimir torque between two uni-axial birefringent dielectric plates has been studied 
more recently in [s*] . The lateral Casimir force between corrugated parallel plates was calculated perturbatively in f3] ■ 
Lately a non-contact rack and pinion arrangement has been proposed in [5| and discussed in the proximity force 
' approximation (PFA) limit. This proposal has been generalized to the design of a non-contact gear consisting of two 
^ ' corrugated concentric cylinders in [6| to discuss possible experimental arrangements. The perturbative calculation for 
the Casimir torque between two concentric corrugated cylinders, in the spirit of 4], has not yet been carried out. We 
^ ^ achieve this task to the leading order for the scalar case here. The next-to-leading-order calculation for the corrugated 
. cylinders, in the spirit of the corresponding case of corrugated plates 0, is in progress. 

o 

oo 
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A. Statement of the problem 

We consider two concentric, semi-transparent, corrugated cylinders described by the potentials, 

V,{r,e)^K5{r-a,~hde)), (1) 

yZ ^ where i = 1, 2 are labels that identify the individual cylinders, and we shall have a — a2 — ai > 0. The functions hiiO) 
describe the corrugations associated with the cylinders. We define the function 

a{e) = a + h2{e)-hi{e), (2) 

which measures the relative corrugations between the cylinders. We shall define the corrugations hi{0) such that the 
mean of the relative corrugations evaluates to a. 



2 cStt 



de[h2{e)~hi{e)\^Q. (3) 
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Let E be the total Casimir energy associated with the two concentric corrugated cylinders, including the divergent 
contributions associated with the single cylinders. In general this energy changes if we rotate one of the cylinders 
with respect to the other and thus leads to a torque between the cylinders. This shift can be described by an angular 
rotation, 9o, of the corrugations on the inner cylinder as hi{9 + 0q) and the corresponding torque will be 



r 



dE 



(4) 



We observe that due to rotational symmetry there will be no torque between two uncorrugated cylinders. Further, 
using the same argument there will be no torque when only one of the cylinders has corrugations on it. In light of 
these observations it is helpful to write the total Casimir energy for two corrugated concentric cylinders in the form 



E = E'-"^ + Ei+ E2 + EuiOo), 



(5) 



where, i?'"' is the energy of the configuration when the corrugations are absent on both the cylinders, Ei is the 
additional contribution to the Casimir energy when one of the cylinders is uncorrugated, and £^12 is the contribution 
to the energy which is present only when both cylinders are corrugated. Thus, E12 is the interaction energy due to the 
presence of corrugation, and only this term in eq. ([5]) contributes to the torque between the cylinders. The physical 
quantities associated with the uncorrugated cylinders thus act as a background, and a reference, and we shall find it 
convenient to denote them by the superscript (0) to mean zeroth order. The potential for the background is 



y/°'(r) =A,5(r-a,), 



(6) 



which has no angular dependence. 

Using the multiple scattering formalism, see 
case of corrugated plates in we can write 



and references in Q, which has been further discussed for the 



AE = E1+E2 + E12, (7) 

where in accordance with the notation given in f^, AX represents the deviation of the physical quantity, X , from the 
background. The term in eq. ([7]) which contributes to the Casimir torque, to the leading order, is given by [7| 



12 2r 



(1) 



where the Green's function associated with the background satisfies the differential equation, 



(8) 



(9) 



In the leading order the corrugations are treated as small perturbations, ^ a < ai, which lets us approximate 

the potentials by 



AV.{r,0) « Vl'\r,0) = h,{0)^vl'\r), 
where we have used the superscript (1) to represent the first order perturbation in the quantity. 



(10) 



II. GREEN'S FUNCTION 

We observe that the evaluation of the interaction energy, £'12 , to the leading order, involves solving for the Green's 
function for the configuration involving the background alone. This amounts to solving the differential equation in 
eq. whose solution can be written as 



(11) 



where = ~ uP' . The reduced Green's function, gin (r^ r'; k), satisfies the equation 

- --^r- - iP - \i5[r - ai) - \25{r - a2) g^^(r,r';K) = 

r or or 



5{r - r') 



(12) 
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JO) 



(r,r';K) = < 



The solution for gm\r,r'; k) in the above equation is expressed in terms of the modified Bessel functions as 
Im{K.r<)K,n{nry) - \ [AiaiA2a2 K1K2 (hKi - I1K2) 

+\iaiKl \2a2Kl] Im{nr)Ira{Kr'), if (r, r') < ai < a2, 
Im{Kr<)K,n{nry) + s [ - A2a2 (l + Aiai hKi)] I,n{nr)I„i{Kr') 
+ \[- Aifli II (1 + Aafla hK2)] K„,{Kr)K^{Kr') 
+ ^ [AiaiA2a2 IiK^] I,n{nr)K„^{Kr') 

+ ^[AiaiA2a2/i-ftr|] KmiKr)Im{Kr'), if ai < {r,r') < 02, 

I,n{Kr^)K,n{Kry) - i [AiaiA2a2 {hKi - I1K2) 

+Aiai II + A2a2 /|] K,n{Kr)Km{Kr'), if ai < 02 < (r,r'), 
i [1 + A2a2 /2i^2] Im{'ir)K^{Kr') - i[A2a2i^|] /™(Kr)/„(Kr'), 

3- I„r{Kr)Km{Kr'), 

^^[l + XiaihKi] I^{Kr)K„^{Kr') - ^[Xiai if] K^{nr)K„,{nr'), 
where we have used the abbreviation 

A = 1 + Aioi hKi + \2a2 I2K2 + AiaiA2a2 hK2 {hKi - I1K2) 



(13) 



if r < fli < r' < a2, 



if 
if 



r < fli < 02 < r 



fli < r < 02 < r 



(14) 



We have used the notation /1.2 = Im{i^0'i,2) and Ki_2 = -R'm('fii.2)- In the other regions not quoted above the sohition 
is determined by using the reciprocal symmetry, gm\r,r'] k) = gm\r',r;K), in the Green's function. 
Using the above expressions we can evaluate 



(15) 



The relevant first derivatives are evaluated using the averaging prescription described in [7,], which is not necessary 
in either the Dirichlet or weak limit, as 



|:5^°H^,r';«) 



r— ai ,r'— a2 



r—a2 ,r'—ai 



K 

A 

K 

A 



I[K2 + ^hK2 

Ik. 

IlK'2-^hK2 



(16a) 
(16b) 



where we have used a prime to denote the derivative of the modified Bessel function with respect to the argument. The 
derivatives acting on the second variable in the Green's can be deduced using the symmetry of the Green's function. 
The relevant second derivatives are evaluated to be 



^ ^ (0)/ I \ 



r—ai ,r'— a2 



K 

a" 



2k 2k, 



(17) 



The above evaluations used the Wronskian, [I„i{x)K'^^{x) — I^^{x)K,n{x)] = —1/x, satisfied by the modified Bessel 
functions. 



III. INTERACTION ENERGY 

In terms of the reduced Green's function defined in eq. (I13p we can write the interaction energy, to the leading 
order, in eq. ([8|) as 



m— — 00 71^ — 00 



(2^) 



(^l)m-n(^2)n-m imriJ (1^) 
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where {hi)m are the Fourier transforms of the functions hi{9), which describe the corrugations on the cyhnders, and 
are defined as 

ih)m= / d9e-'"''^h,i0). (19) 
Jo 

The matrix (or the kernel) Lmli in eq. (jlSp is expressed in the form 

^mn = -^J^ KdKl^l{ai,a2;K), (20) 



where — ^ lo"^ = /c^ + after switching to imaginary frequencies by a Euchdean rotation, lo — > iQ. The related 

(21) 



f 2 

matrix Imnidi, 0-2', n) is expressed as derivatives of the reduced Green's function in the form 



d d [ 

^mil(ai , 02 ; k) = Ai A2 ^ ^ [ r f g^^^ (r, f ; k) 5^°^ (f , r; k) 



r— ai ,r— a2 

The reciprocal symmetry in the Green's function leads to the following symmetry in the above kernel: 

4?n (oi , 02 ; k) = 4m (02 , Ol ; k) . (22) 
The expression in eq. ([2T|) is evaluated using eqs. (fTS]) . (fT6|) . and (fT7| : 



4m(ai:a2;'«) = 



hK2 iiK2 + «:ai /ii^2 (/(i^2 + ^ /ii^2) + «ai (/;i^2 + ^ /ii^2) hK2 
+ na2 hK2 (lik'2 - ^ hK2) + ^tfla [hK!^ - ^ /ii^2) /ii^2 
+ Kaina2 (hK'2 - ^ /1X2) (/( if2 + ^ /li^; 
+Kai^ia2(l[K2 + ^ /iif2) (/ii^2 - ^ Ai^2 



naina2(l[K'2 + ^ hK'2 - ^ /(i^2 - /ii^2) /ii^2 

V 2k 2k 2k 2k / 

KaiKa2/lK2f/(A'2 + ^hK'2 - ^ i'iK2 - ^^hK2 
\ 2k 2k 2k 2k 



(23) 



where we have used the notation in which the modified Bessel function with a tilde on it is of order n and that without 
a tilde is of order m. The tilde on A means that we use the corresponding modified Bessel functions in eq. HH). 

A. Dirichlet Limit 

For the case of the Dirichlet limit (aAi,2 ^ 1) the expression in eq. takes the relatively simple form 

j{2)D^ .^11 1 11 1 (24) 

' ' aia2 [I2K1 - I1K2] [12X1 - iiK2] aia2 Dm{a; kR) Dn{a; kR) 

where the superscript D stands for Dirichlet, and for convenience we have introduced the function 

Dm{a;x) = Im{x[l + a])K,n{x[l - a]) - - a])Kjn{x[l + a]), (25) 

and two related variables: R = (ai + a2)/2, which is the mean radius of the two cylinders under consideration, and 
a — a/2_R, which by definition is less than unity. We note that 



aia2 (1 — a^) 



(26) 



Two cylinders with very large radius, such that ai 00 with a kept fixed, will simulate a parallel plate in the 
region of small variations in the angle {0 — > 0). This corresponds to TO,n — > 00, such that m/R is finite and a — > 0. 
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These limits are compatible with the leading uniform asymptotic approximants to the modified Bessel functions for 
large orders, see for eg. [To| . 



/„(mz)~ ^/^e'"''^^) and K,n{mz) ^ ^ ^e'^^^^^'l m^oo, (27) 



where 



i=-== and ,y(z)^yrT^ + ln— — ==. (28) 



Using the above asymptotic behaviors and neglecting terms of order a we can deduce, for large order m, 

4 a 1 1 



(1 — Q^) Dm{cx; kR) Dn{a; kR) sinh HmCi sinh K„a ' 



(29) 



where we have denoted = + {m/R)"^ and = + [n/R)'^. Using the above limiting form in eq. ([24]) . and 
after interpreting m/ R k as the Fourier transform of the coordinate containing the corrugations on the plates, we 
reproduce the expression derived for the corrugated plates in Q ■ 

Using eq. ((24| in eq. (|20p the i-matrix in the Dirichlet limit takes the form 

aMTryo aia2 £>„(«; Ki?) £>„(«; kJ?)' ^ ' 

which also leads to the corresponding result for the corrugated plates. 



B. Weak coupling Limit 

For the case of weak coupling limit (aAi_2 ^ 1) the expression in eq. ((23)) takes the form 

d d y 1 
-fm«^(ai'°^2;K) = AiA2^ai— 02 lm{K-ai)Km{,na2)In{nai)Kn{Ka2) , (31) 

where W stands for weak coupling limit. Using the above expression for /-matrix in eq. (|20p we can write the i-matrix 
in the weak coupling limit as 

r(2)w _ ^1^2 d d I /ai\ 
47r aai 002 02 Va2 / 

where we have introduced the function 

/>oo 

FmniP) = / X dx K„,ix)K„{x)IrniPx)IniPx). (33) 



Taking the uniform asymptotic approximants of the modified Bessel functions, see eq. (j27p . in the above expression, 
we can reproduce the corresponding result for the corrugated plates in [7] . 

It is possible to convert the above integral into a single sum using the technique described in Q . We replace those 
modified Bessel functions which are well defined at the origin, Im, with their power series expansions, then perform 
the integral using 



CO 







xdxx"^+-+''^KMK,,ix) ^ ^ ; v/r^- wtl ! ->0,n>0,fc>0, (34) 

2 k\(k + m)\(k + ny.(k + m + n)\ 



which leaves eq. (j33p in terms of two sums. One of the sums can be performed after regrouping the terms, and leads 
to 

,2fc 1 kl{k + my.{k + ny.{k + m + 71)! 



(to + n + 2fc + 1)! k'l{k' + to)!(A: - fc')!(fc - fc' + n)! 

1 o2k-\-7n-\-n 

-V— -, TO>0,n>0. (35) 

2^ (2A: + TO + n + l) - ' - ^ ' 

k—Q 



FIG. 1: Non-contact gears: Concentric corrugated cylinders with the same corrugation frequency, v — 15, on each cylinder. Oo 
is the angular shift between the gears. 



Substituting the above expression into eq. (j32p . and taking the derivatives with respect to ai and 02, we immediately 
perform the sum, leading to 



A1A2 1 d r/3™+"+i1 A1A2 1 2 d 

— ■ a 



1-/32 



167r dc 



1 / I - 



a\l + a 



(1-a^) 



TO > 0, n > 0, 



(36) 



where we have denoted (3 — 01/02 for convenience. Interpreting m/R ki and n/R — > k2, and taking the limit 
TO 00, n — > 00, while keeping fci.2 fixed, we obtain the expression for the L-kernel for corrugated plates in the weak 
limit J2- 



IV. SINUSOIDAL CORRUGATIONS 

For the particular case of sinusoidal corrugations, as described in figure [H we will have 

hi{e) = /iisin[i/(6' + 6lo)], (37a) 
h2{0) = /i2sin[i/6i], (37b) 

where /ii 2 are the corrugation amplitudes and v is the frequency associated with the corrugations. Necessarily, v 
must be a positive integer. The Fourier transforms, {hi)m,, corresponding to the above corrugations are 

{hi)m = hi'^\e''''Hra,, - e-^'^''"^™,^ J . (38) 
Zl L J 

In general the corrugation frequencies of the two cylinders can be different. However, we note that, to the leading 
order, the interaction energy gets contributions only when both cylinders have the same frequency. 

Using the above expression in eq. I|18p and using the symmetry property of Imn in eq. ([H]) , which further lets us 
deduce Lmn = Lnm, we can write 

^(2) h h °° h h °" C°° 

-J^ ^ cos{i^Sq) L^^„,j^^ = -cos{u9o)-^ ^ / KdKZ^J^^„_^^(ai, 02; k) (39) 

m— — 00 m— — oo ^ 

where the kernel has been explicitly evaluated in eq. ([23]) . 



Dirichlet Limit 



In the Dirichlet limit the interaction energy in eq. (|39p can be expressed in the form 



=-s(.^?o);^^^5?^"(a) (40) 



2t:RL, ' " 240 a3 a a 



2 4 6 8 10 to 



FIG. 2: Dirichlet limit: Plots of B^'d^;) versus to, for to < 1v and fixed v. The dashed curve is the corresponding plot for 
corrugated plates which is approached by the corrugated cylinders for larger values of v. 



where we have divided by a factor of 27ri?, which is the mean circumference in the direction of corrugations. We have 
also defined a suitable function Bi, (a) to make it convenient to compare our results with those obtained in the PFA 
limit and with those for corrugated plates in the appropriate limits. We define 



8a" 



: dx 



1 



1 



(1 - q;2) Dm{a; x) Dm+u{a; x) ' 



(41) 



Using eq. ([^5]) it is straightforward to verify the corrugated plate limit of the above expression. The function B^'^^ {a) 
has been plotted with respect to to = lav in figure [2] The redefined parameter helps us compare our results with 
the corrugated plates. We note that for larger values of v the plots approach the curve for the corrugated plates very 

quickly. We note that B')f'^ [\) — 0, because for a = 1 the radius of the inner cylinder approaches zero. Thus, it is 
pointless to consider the regime a > 1. We also note that B^^^{0) — 1, which then implies the PFA limit. 

The Casimir torque per unit area, for the Dirichlet case, can thus be evaluated, using eq. ((40)) in eq. (|4|), to be 



-j-{2)D 



240 a- 



(42) 



B. Weak coupling limit 



The evaluation for the L-matrix in the weak limit in cq. (|36p is valid for positive indices only. Therefore, we begin 
by rewriting the expression for the interaction energy in eq. (|39p in the form 



E'. 



{2)W 



cos(i/0o) 



hih2 



„ j-{2)W 



j(2)W 



(43) 



where the finite sum is interpreted as zero when u = \. We have used the symmetry property mentioned before 



eq. (P^ . and further used Lm|^ = L^^l^\i — Lm'"n — ^-m.-n^ which can be deduced from eq (PT|) using the Bessel 

into 



-(2)W 



-(2)W 



-(2)W 



function property /_m(a;) = Im{x) and K^m{x) = Km{x). After substituting the L- matrix, derived in eq. ((36 
the above equation we can immediately perform the sums to yield 



E. 



[2)W 



2ttRL^ 



cos (j^ 6*0) 



A1A2 



327r^ a a a 



(44) 



where we have defined the function 



' 2 da 



1 /I 



a^ 



a 



2 ll + a 



{I - a^){l + 2ay + a^) 



(45) 
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FIG. 3: Weak coupling limit: Plots of bI^^^ {^) versus to, for to < 1v and fixed v. The dashed curve is the corresponding 
plot for corrugated plates which is approached by the corrugated cylinders for larger values of v. 



The Casimir torque per unit area, for the weak couphng Hmit, can thus be evaluated, using eq. 

j/sin(j/6'o) —T, (a). 



2-kRL, 



327r2 a a a 



in eq. ([3]), to be 

(46) 



We note that B, 



■2)W , 



Dirichlet case, we note that B^^'^ [1) 



(0) — 1. This verifies that the above result satisfies the proximity force theorem. As in the 

0, because for a = 1 the radius of the inner cylinder approaches zero. The 
above result should also yield the result for corrugated parallel plates in the limit u —> oo, oi _2 — > oo, i? cxd, such that 
a and u/ R is finite. Also, in this limit a ^ 0. Recalling the corrugated plates parameter [3], k^a va/ R = 2i>a = to, 
we note that the limit to corrugated plates is achieved by taking the limit — > oo with to kept fixed. To this end we 
rewrite eq. (|45|l in terms of to as 



2 dto 



1 



2v 



1 



in which the 



oo limit can be immediately taken to yield 
to 



tl^ 

2 dto 



^{l + to)e- 



to 


t^ 








2 


dtl 





(47) 



(48) 



which matches the result for the corrugated plates exactly. The function B'if''^ {a) has been plotted with respect to 
to in figure [31 As in the case of Dirichlet case, we note that for larger values of v the plots approach the plot for 
corrugated plates very quickly. 



V. CONCLUSION AND FUTURE DIRECTIONS 



We have evaluated the Casimir torque between two concentric corrugated cylinders perturbatively in the corrugation 
amplitude for the scalar case. We have derived explicit expressions for the case when the cylinders have sinusoidal 
corrugations. Nonzero contributions in the leading order requires the corrugation frequencies on the two cylinders 
to be identical. Our results for the Casimir torque reproduce the results for the lateral force on corrugated parallel 
plates in the limit of large radii and small corrugation wavelengths. 

Extension of the calculation to the next-to-leading order in the spirit of ,7] is in progress. We recall that it was 
possible to evaluate the lateral force between corrugated parallel plates, in the weak limit, exactly in terms of a single 
integral 0. The corresponding exact result for the Casimir torque between corrugated cylinders, in the weak limit, 
is being sought. Generalization to the electromagnetic case is most important and we intend to attempt it next. 
Generalization of the calculation to the case of a rack and pinion arrangement is another possible extension. We hope 
that these calculations will have application in the design of practical nanomechanical devices. 
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